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Abstract 

This paper is a direct continuation of [|l|] where we begun the study of the integrable 
structures in Conformal Field Theory. We show here how to construct the operators 
Q±(A) which act in highest weight Virasoro module and commute for different values of 
the parameter A. These operators appear to be the CFT analogs of the Q - matrix of 
Baxter 0, in particular they satisfy famous Baxter's T — Q equation. We also show that 
under natural assumptions about analytic properties of the operators Q(A) as the functions 
of A the Baxter's relation allows one to derive the nonlinear integral equations of Destri-de 
Vega (DDV) for the eigenvalues of the Q-operators. We then use the DDV equation 
to obtain the asymptotic expansions of the Q - operators at large A; it is remarkable that 
unlike the expansions of the T operators of 0], the asymptotic series for Q(A) contains the 
"dual" nonlocal Integrals of Motion along with the local ones. We also discuss an intriguing 
relation between the vacuum eigenvalues of the Q - operators and the stationary transport 
properties in boundary sine-Gordon model. On this basis we propose a number of new 
exact results about finite voltage charge transport through the point contact in quantum 
Hall system. 



1 



1. Introduction 



Existence of an infinite set of mutually commuting local Integrals of Motion (IM) is 
the characteristic feature of an integrable quantum field theory (IQFT). Therefore simul- 
taneous diagonalization of these local IM is the fundamental problem of IQFT. In the 
case of infinite-size system this problem reduces to finding mass spectrum and factorizable 
S-matrix associated with IQFT; much progress in this direction has been made during the 
last two decades (see e.g. for a review). On the other hand, for a finite-size system 
(say, with the spatial coordinate compactified on a circle of circumference R) this problem 
becomes highly nontrivial and so far its solution is known to a very limited extent. Most 
important progress here has been made with the help of so called Thermodynamic Bethe 
Ansatz (TBA) approach 0,0. TEA allows one to find the eigenvalues associated with the 
ground state of the system (in particular the ground-state energy) in terms of solutions of 
nonlinear integral equation (TBA equation). However it is not clear how the combination 
of thermodynamic and relativistic ideas which is used in traditional derivation of the TBA 
equation can be extended to include the excited states. 

The above diagonalization problem is very similar to that treated in solvable lattice 
models. In the lattice theory very powerful algebraic and analytic methods of diagonal- 
ization of the Baxter's families of commuting transfer-matrices are known [0, these 
methods are further developed in Quantum Inverse Scattering Method (QISM) Of 
course many IQFT can be obtained by taking continuous limits of solvable lattice models 
and the method based on commuting transfer-matrices can be used to solve these QFT. 
This is essentially the way how IQFT are treated in the QISM. However, for many IQFT 
(notably, for most of IQFT defined as perturbed CFT |T^]) the associated solvable lattice 
models are not known. Besides, it seems to be conceptually important to develop the 
above methods directly in continuous QFT, in particular, to find continuous QFT versions 
of the Baxter's commuting transfer-matrices. 

This problem was addressed in our recent paper where we concentrated attention 
on the case of Conformal Field Theory (CFT), more specifically on c < 1 CFT. We should 
stress here that although the structure of the space of states and the energy spectrum in 
CFT are relatively well understood the diagonalization of the full set of the local IM remains 
very nontrivial open problem. In we have constructed an infinite set of operator valued 
functions Tj(A), where j = 0, and A is a complex variable. These operators 



(we will exhibit their explicit form in Sect. 2) act invariant ly in irreducible highest weight 
Virasoro module Va and they commute between themselves for any values of A, i.e. 

T,(A): Va^Va, 

(1.1) 

[T,(A), T,'(A')] = 0. 

The operators Tj(A) are defined in terms of certain monodromy matrices associated with 
2j + 1 dimensional representations of quantum algebra Uq{sl2) where 

q = e'""^^ . (1.2) 

and j3 is related to the Virasoro central charge as 

c = 13 - 6 + . (1.3) 

Evidently, the operators Tj(A) are CFT versions of the commuting transfer-matrices of 
the Baxter's lattice theory. We will still call these operators "transfer-matrices" although 
the original meaning of this term apparently is lost. As we have shown in [0, in CFT 
the operators Tj(A) enjoy particularly simple analytic properties, namely they are entire 
functions of A^ with an essential singularity at A^ = oo and their asymptotic behavior 
near this point is described in terms of the local IM. Therefore the operators Tj(A) can 
be thought of as the generating functions for the local IM since all the information about 
their eigenvalues is contained in the eigenvalues of Tj(A). The operators Tj(A) are shown 
to obey the "fusion relations" which for any rational value of in ( p..2|) provide a finite 
system of functional equations for the eigenvalues of these operators. For the ground-state 
eigenvalue (in CFT it corresponds to a primary state) these functional equations turn 
out to be equivalent to the TBA equations; in general case they provide modified TEA 
equations suitable for the excited states. Interesting but somewhat inconvenient feature of 
this approach is that the resulting TBA equations depend on c in a very irregular manner 
(they depend on the arithmetic properties of the rational number 0^) whereas the resulting 
eigenvalues of Tj (A) are expected to be smooth functions of c. 

Another powerful method known in the lattice theory is based on the so-called Q 
-operator. This method was introduced by Baxter in his original study of 8-vertex model 
0. One of its advantages is that it is not limited to the cases when g is a root of unity. 
In this paper (which is a sequel to 0) we introduce the analog of Q-operator directly in 
CFT and study its properties. The Q-operators (in fact we will define two Q-operators, 
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Q-I-(A)) are defined again as the traces of certain monodromy matrices, tliis time associated 
with infinite- dimensional representations of so called "g-oscillator algebra" . The operators 
Q±(A) obey the Baxter's functional relation 

T(A)Q(A) = Q(gA) + Q(g-U) , (1.4) 

where T(A) = Ti (A). This construction is presented in Sect. 2 where also the most impor- 
tant properties of the Q-operators are discussed. 

In the lattice theory the Baxter's relation (|1.4| ) is known to be a powerful tool for 
finding the eigenvalues of the transfer- matrices , the knowledge about analytic properties 
of the Q-operator being a key ingredient in this approach. Our construction of the Q- 
operators as the traces of the monodromy matrices makes it natural to assume that they 
enjoy very simple analytic properties similar to those of the operators Tj(A): up to overall 
power-like factors they are entire functions of A^ with the following asymptotic at A^ —oo 

log Q±(A) ~ M (-A^)^^ , (1.5) 

where M is a constant (which will be actually calculated in Sect. 3). Using these properties 
we show that the eigenvalues of the Q-operators satisfy the Destry-de Vega (DDV) equation 
1^. This is done in Sect. 3. The DDV equation can be solved exactly in the limit A 
+00, where A is the Virasoro highest weight in (|1 . Ij ) , and we analyze in the Sect. 3 the 
vacuum eigenvalues of the operators Q(A) in this limit. In Sect. 4 we further study the 
properties of the Q - operators and formulate our basic conjectures about their analytic 
characteristics. The exact asymptotic expansions of the Q and T operators at A^ ^ oo 
are proposed here. We observe that unlike the asymptotic expansion of T(A) the large 
A^ expansion of Q(A) contains both local and nonlocal IM and that the operators Q(A) 
obey remarkable duality relation with respect to the substitution —>■ Although the 
results of this section have somewhat conjectural status we support them by explicit study 
of the eigenvalues of the Q operators at the "free fermion point" = 1/2. In Sect. 5 we 
discuss the relation of the Q operators to the characteristics of stationary non-equilibrium 
states in so called boundary sine-Gordon model [|ll|] (see also [|12|); these states attracted 
lately much attention in relation to the finite-voltage current through the point contact 



in a quantum Hall system [|T^, [0, [jT5|, [0. Possible directions of further studies are 
discussed in Sect. 6. 
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2. The Q-operators 

In this section we will introduce the operators Q±(A) which satisfy ( p..4|) . We start 
with a brief review of the definitions and results of [Q. 

Let ip{u) be a free chiral Bose field, i.e. the operator-valued function 



n 



(2.1) 



Here P, Q and an,n = ±1, ±2, . are operators which satisfy the commutation relations of 
the Heisenberg algebra 



Z 2 ^2 



(2.2) 



with real 13. The variable u is interpreted as a complex coordinate on 2D cylinder of a 
circumference 27r. As follows from (|2.1|) the field ip{u) is a quasi-periodic function of w, 
i.e. 

(/?(w 4- 27r) = (/?(w) + 27rzP . (2.3) 

Let jFp be the Fock space, i.e. the space generated by a free action of the operators a„ 
with n < on the vacuum vector | p) which satisfies 



0"n \ p) = ^ , for n > ; 
P \p) = p \ p) ■ 



(2.4) 



The composite field 



(2.5) 



is called the energy-momentum tensor; it is periodic function of u and its Fourier modes 



du 
2^ 



(2.6) 



generate the Virasoro algebra with the central charge (|1.3|) [|T7| , [|T8| . It is well known 
that for generic values of the parameters (3 and p the Fock space Tp realizes an irreducible 
highest weight Virasoro module Va with the highest weight A related to p as 



P\^ c- 1 
(3) 24 



(2.7) 



For particular values of these parameters, when null- vectors appear in jFp, Va is obtained 
from Tr> by factoring out all the invariant subspaces. The space 



-r- _ moo 



D 'p 



admits the action of the exponential fields 

V±{u) =: e±2^(") : . (2.9) 
Also, let E, F and if be a canonical generating elements of the algebra Uq(^sl{2)) i.e. 

[if, E] = 2E, [H, F] = -2F, [E, F] = , (2.10) 

where q is given by (|1.2|) . Let j be a non-negative integer or half-integer number. We 
denote ttj an irreducible 2j + 1 dimensional matrix representation of t/g(s/(2)) so that 
Ej = 7:j{E), Fj = T^j{F) and Hj = T^j{H) are (2j -|- 1) x (2j + 1) matrices which satisfy 
the relations ( |2.1C1| ). 

The "transfer-matrices" Tj(A), j = 0, |, 1, |, ... are defined as 



oo 



T,(A) =tr, 

X) 



2-K 



Pexp(A / Kj{u)du) 



n=0 



Kj{ui)Kj{u2)---Kj{Un) duidu2---dur 



27r>lti >112 > . ..>Un >0 



(2.11) 



Here 

Kj{u) = V-iu) q^Ej + V+{u) q'^Fj 
and V denotes the operator ordering along the integration path. Obviously, 

To(A) = I , 



(2.12) 



where I is the identity operator. Although the exponentials in ( p.ll|) act from one com- 
ponent of the sum ( |2.8|) to another it is easy to see that the operators (|2.11 ) invariantly 
act in the Fock space J^p. Another obvious property of Tj is that these operators are in 
fact the functions of A^ as the traces of the odd-order terms in ( p.ll| ) vanish. 

As is explained in [|I| the operators ( |2.11| ) form the commuting family, i.e. they satisfy 
the relations ( |1 . 1|) . They also obey the "fusion relations" 



T{X)T,{q^+H) = T^ i (g^+U) + T^+i (?^A) . 



^2.13) 
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Together with (|2.12 ) this relations aUow one to express recurrently any of the "higher-spin" 
operators Tj(A) with j = 1, |, 2, ... in terms of the basic one 



T(A) = Ti(A) . 



(2.14) 



In the case j = ^ it is easy to evaluate the traces in ( |2.11j ) and obtain the power-series 
expansion 

oo 

T(A) = 2 cos(27rP) + ^^"^ , (2.15) 



n=l 



where the coefficients define an infinite set of basic nonlocal IM 



r / ( e'^"^ V.{ur)V+{u2)V.M...V+{u2n) + 

'27r>ltl>U2>...>U2n^0 



-2i7rP 



V+{ui)V-{u2)V+{uz).V-{u2n) ) dUidU2...dU2n ■ 



It follows from (|1 . 1|) that these nonlocal IM commute among themselves 



[G^, Grn\ — . 



(2.16) 



The operator-product expansion 

V+{u)V-{u') = {u-u')-'^l^\l + 0{u-u')) , u-u' (2.17) 



shows that the expressions ( |2.11| ) and ( pj.l6D can be taken literally only if 

< < L 



^2.18) 



for otherwise the integrals in (|2.16|) diverge. In what follows we will call the region ( |2.18| ) 
the Semi-classical Domain (SD). In fact one can define the operators G^ outside the SD 
by analytic continuation in /3^. A convenient way to do that is to transform the ordered 
integrals in ( |2.16| ) to contour integrals. For example Gi can be written as 

Gi= (q2_^-2^-1 1 ^^^-2.^P _^-1^2.^P^^^ 

Jo Jo I (2.19) 

V-(ui+iO)V+{u2-iO) + {qe^'''^ -q-^e-'^'''^) V+(ui + iO)V-(u2 - iO) | . 



Obviously, ( |2.19| ) does not contain divergent integrals and for in SD coincides with 
the ordered integral in ( p.l6| ). Similar representation exists for higher G^, and thus the 



operator T(A) can be defined outside SD through (|2.15|) . Of course the operators G„ and 



T(A) obtained this way exhibit singularities at = 



277,-1 

I3l = — , n = 1,2,3,..., (2.20) 



where the integrals ( |2.16|) develop logarithmic divergences. In order to define the operators 



Tj(A) at the singular points ( p.20|) some renormalization is needed but we do not discuss 
it here (see however our analysis of the case (3^ = 1/2 in Sect. 4). 

As is mentioned in the Introduction, the power series ( p.l5|) defines T(A) as an entire 
function of A^ with an essential singularity at A^ — oo. Its asymptotic expansion near this 
essential singularity can be expressed in terms of local IM as 



1-2™ 



logT(A) ~m AT^ l2n-i , (2.21 



n=l 



where l2n-i is the basic set of commutative local IM as defined in [jT|; the operators l2n-i 
can be written as the integrals 



•-271-1 







'''^T2.(«), (2.22) 



where the local densities T2n{u) are particular normal ordered polynomials of du^iu),..., 
d'^ipiu) (see [|l| for details and for our convention about the normalization of l2n-i)- The 
numerical coefficients m and C„ in the expansion (|2.21| ) will be calculated exactly in Sect. 3. 
As follows from ( |1 . 1|) all these local IM commute with the nonlocal IM G^. 

As is known [^, [^, the local IM l2n-i defined in ^ do not change under the 



substitution /S"^ ^ /3 ^, if we simultaneously make the replacement (p{u) — > (3 ^^{u) 

l2n-i{/3',¥^(w)} = l2n-i{/3"',/3"V(«)} . (2.23) 

Evidently the nonlocal IM do change and so there exists an infinite set of "dual" nonlocal 
IM Gn which are obtained from ( |2.16| ) by just this substitution, i.e. 

G„ = / ( f'' U-{ui)U+{u2)U-{us)...U+{u2n) + 

J2Tr>ui>U2>...>U2„>0 (2.24) 

g"^^ U+{ui)U-{u2)U+{us)...U-{u2n) ) dUidU2...dU2n , 

where 

U±{u) =: e 13^^^ ' : 
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and 



q = e^ . (2.25) 



Of course if < 2 the analytic continuation described above is needed to define G^- How- 
ever it is possible to check that the operators thus defined commute among themselves 
and commute with all the nonlocal IM G^ and local IM l2m-i 

[G„, G^] = [G^, G^] = [G^, l2m-i] = . (2.26) 

Now we are in position to define the Q-operators and to describe their basic properties. 
Let S+,S- and Ti. be operators which satisfy the commutation relations of so called "q- 
oscillator algebra", 

qS+£. - q-^£-£+ = , [n,£±] = ±2£± (2.27) 

and let p be any representation of this algebra such that the trace 

Z{p) = trp[ e^^^P^ ] (2.28) 

exists for complex p belonging to the upper half plane, p > 0. Then one can define 
two operators 

A±(A) = Z-i(±P)tr4e±2-^'^Pexp(A / du{V-{u)q^'^ £± + V+{u)q'^'t £^)) ] , (2.29) 

Jo 

where again the symbol V denotes the operator ordering along the integration domain. As 
we are interested in the action of these operators in J-'p the operator P in ( |2.29[ ) can be 



substituted for its eigenvalue p. Strictly speaking the definition ( |2.29| ) makes sense only if 
Qm p > for A_|_ and if 53m p < for A_. However these operators can be defined for all 
complex p (except for some set of singular points on the real axis) by analytic continuation 
in p. Then it is easy to see that A_ can be obtained from A_|_ by substitution 

P -P, ip{u) -ip{u) . (2.30) 

The operators (|2.29| ) can be written as the power series 



A±(A) = 1 + ^ A^" a2n(ai, . . . , a2n\ ± P) J2n{Tcri, TCr2, . • • , Tcr2n) , (2.31) 



n=l {<Ti = ±i} 

CTlH ho-2„=0 



where 



J2n(0-l, . . . ,CT2n) = / Va^{ui)Va2{u2) ■ ■ ■Va2^{u2n) dUl . . .dU2n 

J 27r >Ul • ■■U2n^0 

and 

a„(ai, . . .,a2n\P) = Z-\P) tr p(^e'^''"^£^,£^, ■ ■ ■ £^,^ . (2.32) 

It is easy to see that the coefficients ( |2.32p are completely determined by the commutation 
relations (|2.27| ) and the cyclic property of the trace and so the operators ( |2.29| ) do not 



depend on the particular choice of the representation p. The operator coefficients in the 
power series ( p.31| ) can be expressed in terms of the basic nonlocal IM ( p.l6| ). It is not 
difficult to calculate first few terms of this series 



A±(A) = 1 - J,. - A 



4sinasin(a ± x) [ 4sin2asin(2a ± x 

G? 1 \6 f G 



16sinasin2asin(a ± x) sin(2a ± x) j |^ 4sin3asin(3a ± x 

Gf — 4 GiG2(sinasin(a ± x) + sin2asin(2a ± x)) 



<! 7— ^77— kr^+ (2-33) 
0(A«) , 



64 sin a sin 2a sin 3a sin(a ± x) sin(2a ± x) sin(3a ± x) 
where 

X = 27rP, a = 7r/32 . 

For further references it is convenient to introduce an alternative set of nonlocal IM defined 
as coefficients in the expansion 

oo 

logA+(A) = -^ y'^^H^, (2.34) 

n=l 

where 

y = r^T{l-P^)\. (2.35) 



These coefficients are, of course, algebraically dependent on those in (|2.16| ) . For example 



(3'^ r(/32) 

= 47rr(l-/32) sin(27rP + 7r/32) ' ^^'^^^ 



Define also a new set of "dual" nonlocal IM 



H4/3^(^(^i)} = H4r^/3-v(«)} , 
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(2.37) 



obtained from H„ by the analytic continuation, as described above. 
The operators Q±(A) are defined as 

Q±(A) = A±2^/^'a±(A) . (2.38) 

Like the operators Tj(A) above the operators ( p.38| ) act in a Fock space 

Q±(A): T^^T^. (2.39) 

The operators Q±(A) exhibit remarkable properties. Here we simply list some of them 
leaving the proofs to the other paper. 

i. The operators Q±(A) commute among themselves and with all the operators Tj (A), 

[Q±(A), Q±(A')] = [Q±(A), T,(A')] = . (2.40) 

ii. The operators Q±(A) satisfy the equation (pT^), i.e. 

T(A)Q±(A) = Q±(qA) + Q±(q-U) . (2.41) 

The equation ( p..4| ) can be thought of as the finite-difference analog of the second order 
differential equation so we expect it to have two linearly independent solutions. As T(A) 
is a single- valued function of A^, i.e. it is periodic function of log A^, the operators Q±(A) 
are just two "Bloch-wave" solutions to the equation ( [1.41 ). The operators Q±(A) satisfy 
the "quantum Wronskian" condition 

Q+(QiA)Q_(g-5A) - Q+(Q-n)Q_(gn) = 2i sin(27rP) . (2.42) 

Hi. The "transfer-matrices" Tj(A) can be expressed in terms of Q±(A) as 

2i sin(27rP) T,(A) = Q+(Q^+5A)Q_(g-^'-U) - Q+(g-^'-5A)Q_(g^+n) . (2.43) 

In view of this equation the operators Q±(A) appear more fundamental then the transfer- 
matrices Tj(A). We will see more support to this idea below. 

Let us just briefiy sketch the derivation of the functional relation of this section. 



The main idea is to consider more general T-operators T^(A) defined as in ( |2.11| ), but 
associated with the infinite dimensional representation Tr"*" of t/q(s/(2)) with arbitrary 
(complex) highest weight 2j. Note that if j takes a no n- negative integer or half-integer 
value then the matrices 7r^(i?), 7Tj~{F) and 7rj~{H) have a block-triangular form with two 
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diagonal blocks, one equivalent to the (2j ' + l)-dimensional representation tTj and the other 
being the highest weight representation 7r^^_j^. In this way we obtain the following simple 
relation 

T+(A) = T,(A) + Ti^_i(A) , J =0,1/2,1,3/2,... . (2.44) 
Next, the operator T^(A) enjoy the following remarkable factorization property 



2zsin(27rP) T+(A) = Q+(?^+^ A) Q-iq~'~^ A) 



(2.45) 



which is proved explicitly by using decomposition properties of the tensor product of two 
representations of the g-oscillator algebra (the latter are also representations of the Borel 
sub-algebra of t/g(s/(2)) with respect to the co-multiplication from C7q(s/(2)). Then the 
functional relations ( |2.42| ), ( p.43|) trivially follow from (|2.44|) and ( p.45|) above, while the 
remaining relations ( |2.40| ), ( |2.41| ) are simple corollaries of these two. 

Although our considerations here were specific to the continuous theory similar results 
hold for Q-matrix of the lattice theory as well. This is quite obvious since the structure of 
the functional equations is determined merely by the decomposition properties of products 
of representations of Uq[sl{2)) associated with the monodromy matrices. The details of 
these calculations will be given in ||23|| . 

The operators Q±(A) take particularly simple form when applied to the space J-'p 
with 2p equals some integer which we denote A^. In these cases the "quantum Wronskian" 
( p.42|) is equal to zero and the solutions A_|_(A) and A_(A) ( |2.38|) coincide . For 2p = N 
all the coefficients (|2.32| ) are readily calculated 



a2n(cri, . . . , o-2n| A/2) = {q - q' 



-2n 



(2.46) 



and the operators (|2.38|) can be written as 



A±(A) 
where 



p=N/2 



E 



n=0 ^ ' 



q'^V 



27r 



V+{u) 



ttA 



2-K 



du 
2^ 



V-{u) 



p=N/2 



, (2.47) 



(2.48) 



sin(7r/32) 

and, as before, the symbol V denotes "^-ordering" , i.e. the operators V±{u) with greater 
u are placed to the left. For p = A/2 the exponentials V±{u) are 27T-periodic functions of 
u and so the integration contours in ( p.47| ) close. The series ( p.47| ) is closely related to 
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the boundary state in so called Boundary Sine-Gordon model (with the bulk mass equal 
zero)[jll|. In particular, its vacuum- vacuum matrix element A^]^'^'^\x) , defined as 

A^iX) \ p) = A^r^\X) \ p) , (2.49) 

coincides for p = N/ 2 with the one-dimensional Coulomb gas partition function 



2tt 



n 

X 



2 0/32 



Jj|4sin(^^i^)sin(^^)f' n|2«-(^^^) 



20" 



(2.50) 

As was shown in [ p^ the series ( p.50| ) defines an entire function of p,'^ with the asymptotic 
behavior logZN{p) ~ const {p^) 2-20^ . In fact, it is easy to show that this result implies 
not only to the vacuum eigenvalue (|2.49| ) but also to any matrix element of the operator 
A. We conclude that for p = N/2 the operators A±(A) are entire functions of and they 
enjoy the asymptotic form 

log A±(A) ~ M (-A^)^^ , A^^-cx). (2.51) 



We should stress again that the relations ( |2.47| ) and ( p.50|) hold only for 2p = N . For 



non-integer 2p they do not hold. Nonetheless we find it natural to assume that the above 
analytic properties of the operators A±(A) as the functions of A^ hold for any p. In the 
following Sections we use the functional equation ( |1.4|) together with this analyticity 
assumption to derive various asymptotic expansions for the operators A±(A). 



3. Destri-de Vega equation 

Now let us turn to the eigenvalue problem for the operators Q±(A) 

Q±(A) |a) = Q^(A) |a) , (3.1) 

where | a) G Tp. In this section we consider only the case of 0^ in SD (|2.18|) . Let us 
concentrate attention on one of the operators ( |2.38| ), say Q(A) = Q_|_(A); the problem for 
Q_(A) can be solved then by the substitution (|2.30| ). 
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As follows from ( 2.41 ) any eigenvalue Q{X) of Q(A) satisfies the Baxter's functional 
equation 

TiX)QiX) = QiqX) + Qiq-'X) , (3.2) 

where T(A) is the corresponding eigenvalue of the operator T(A). Denoting A{X) the 
eigenvalue of the operator A_|_(A) in ( |2.3(j| ) one can rewrite ( P?^ ) as 

T(A)A(A) = e^'^'PAiqX) + e'^'^'P A{q-^ X) . (3.3) 

Although this equation looks as just a relation between two unknown functions in fact it im- 
poses severe restrictions on these functions provided their analytic properties are known0. 
Motivated by the relation to the Coulomb gas partition function ( |2.5U| ) discussed in Sect. 2 
we accept here the following assumptions about analytic properties of the eigenvalues A{X). 
Let (3'^ be restricted to SD and let 53m p = 0. Then 

(i) Analyticity. The functions A{X) and T(A) are entire functions of the complex 
variable A^ i. 

(ii) Location of zeroes. Zeroes of the function A{X) in the A^-plane are either real or 
occur in complex conjugated pairs. For any given eigenvalue A{X) there are only finite 
number of complex or real negative zeroes. Real zeroes accumulate toward -foo in the A^. 
For the vacuum eigenvalues all zeroes are real and if 2p > — they are all positive. 

(iii) Asymptotic behavior. The leading asymptotic behavior of A{X) for large A^ is 
given by ( p.51|) with some constant M. 

If < < ^ an entire function A{X) with the asymptotic behavior (|2.51|) is com- 
by its zeroes A^, 

product 



pletely determined by its zeroes A?, /c = 0, 1, . . . ; it can be represented by a convergent 



.(A).nO-^). (-) 

where the normalization condition 

A{0) = 1 (3.5) 

is taken into account. 

^ In lattice theory the equation ( |3.3| ) (which appears there "decorated" with some non-universal 
factors in the right hand side) leads to the Bethe Ansatz equations which completely determine 
the eigenvalues [Q. 

^ In fact one can prove this property for T(A), see Q. 
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Actually, one important remark about the normalization (|2.38| ), ( p.5|) of the operators 
A must be made here. For given eigenvalue A{X) the positions of its zeroes Xk depend on 
p and for special values of this parameter one of these zeroes, say /cq, can happen to be 
at zero, Afco = 0. For example, for the vacuum eigenvalue A^_^"'^\x) this happens at 2p = 
—(3'^—n] n = 0, 1, 2, as brief inspection of ( p.33|) and (3.21) below shows. Obviously, at 



these values of 2p the normalization of the operators A as in (|2.38| ), ( |3.5| ) is not suitable 
as under this normalization all the terms in the power series A{X) = 1 + Xl^i ^^^"^ 
except the first one would diverge. Of course this formal singularity can be eliminated by 
renormalization of A with appropriate p-dependent factor. As our results below are not 
sensitive to this subtlety, in what follows we use the normalization ( |2.38[ ), ( ^.51 ) (which we 
find technically very convenient) and just assume that p does not take these "dangerous" 
values (i.e. that that all 7^ 0). 
Introduce the function 

^ .--^ . (3.6) 

Setting A^ = A| in ( |3.3| ) and remembering that T(A) has no singularities at finite A^ one 
obtains the Bethe-Ansatz type equations for the positions of the zeroes A| 

a(Afc) = -1 . (3.7) 

This is an infinite set of transcendental equations for the infinite set of unknowns A^. It 
can be transformed into a coupled system of a single non-linear integral equation and a 
finite set of equation ( |3.7D only for those zeroes which do not belong to the positive real 
axis in the A^-plane. This system is known as Destry-de Vega (DDV) equation [0. In the 
Appendix A we show that under the assumptions (i)-(iii) given above the infinite set of 
equations ( |3.7| ) is equivalent to 

i log a{e) = + 2Mcos^ + ij^'log 3(9 - Oa) - 2G * log (l + a{e - iO)) 



aiOa) = -1 , 
where we have used new variables 



(3.8) 



J- + ? 
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and denote zeroes of A{X) lying outside the positive real axis of A^. The star m 
denotes the convolution 



/oo 
de'A{e - e')B{e') (3.10) 
- 00 



and 



Q^m ) -r^^ . ^ sinh(7rKl+0/2) ^ 

=exp<^-z / — s\n{u6) — \ . — — } , 3.11 

' Jq V cosh( 7rz^/2 ) sinh( 7rz^^/2 j ' 

G{(^) = m + 7^ ^elog 3(6) . 

Note that the function S{9) coincide with the soliton-soliton scattering amplitude for the 
Sine-Gordon model [|5|. The equation ( p.8|) for the vacuum eigenvalues (in fact, more 
general equation obtained by replacing in the r.h.s. by 2sinh^) was originally derived 
from the Bethe ansatz equations associated with XXZ lattice model by taking appropriate 
continuous limit which leads to the sine-Gordon field theory 0. Here we obtain (the 
"mass-less" version of) this equation directly in continuous field theory, bypassing any 
direct reference to the lattice theory. 

Given a solution of ( |3.(j| ) the function logA(A) can be calculated as 

logA A =-M dv — — — ^TTT (-A ) , 3.12 

Jc^ cosh7rz//2 smh7rz^4/2 ^ ' 

where the function g{v) is defined as 

/+00 in 
— 53m log (1 + 0(6* - zO)) e-*"V (3.13) 
-oo 



The integration contour Cjj in ( |3.12D goes along the line Q'm u = —1 — t with arbitrary 
small positive e. 

It is instructive to study the vacuum eig envalue A(^"^)(A) of A(A), defined by ( p:i9| ), 
in the limit p -foo. Note that the parameter p enters the DDV equation ( ^.8| ) exactly 
the same way the external gauge field (which couples to the soliton charge) appears in 



the sine-Gordon DDV equation [^. The DDV equation ( p.8| ) is greatly simplified in the 
limit p According to the assumption (ii) above all zeroes A^ in this case are real 

and positive and so one can drop the term «X]a'^°§ "^(^ ~ r.h.s. of the integral 

equation ( ^.8| ). Moreover, in the limit p -|-oo this equation reduces to a linear equation 
of Winer-Hopf type (see Appendix A for some details) 

--f + Mcos^ e'- — delog S{9 - 9') Qmlog (l + a(""^)(^')) = , (3.14) 
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provided one assumes that 



1 + ^ 

^ip) = — 7^ — log -^0 ~ const logp , p +00 



(3.15) 



where Aq is a minimal of the zeroes A|. 



The equation ( p.l4| ) can be solved by the standard technique As a result one 

obtains 



,(vac)r.\ iP r(l-izy(l+0/2) 



iSv f \2 \-«!^(l+C)/2 



4z/ r(3/2 - iiy/2) T{1 - iui/2) 
where g^'"°-'^\i>) is related to a'^""^^(6') through ( |3.13|) 

,2 N (1+0/2 



(3.16) 



(3.17) 



and 



25 =(1 + log(l+0-eiog^ 



Note that (|3.17| ) supports the assumption ( p.l5|) . The equation (|3.12| ) gives 
logA("-)(A)|p^+oo- 

^ r(i - + 0/2) r(i + zz.e/2) r(-l + iv/2) e'- ( - ^ 



p 



A2 x^!-(1+C)/2 

2 




(3.18) 



where integration contour C^, go along the line v = — 1 — e. The integrand a"(z^) 

in ( |3.18| ) has the simple asymptotic behavior in whole complex z^-plane except positive 
imaginary axis 



P 



sinh(7rj^(l+0/2) V Ag 



A2 n *t'(l+0/2 



— >■ oo 



and we can work out the integral ( |3.18|) for | A^ | < | Aq | by closing the integration contour 
Cjy at the infinity in the lower half plane u < —1 and calculating the residues of the 
poles at z/ = — 2m/(l +0? n = 1,2, ... . The result can be written in the form 



logA(-^)(A)|p^+oo = - E ^n"'^ 



p — >00 ; 



(3.19) 



n=l 
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where y is given by (|2.35| ) and the above residue calculations give the coefficients 
Hn"''^^\p^oo explicitly. According to (|2.34| ) this coefficients coincide with the large p asymp- 
totic of the vacuum eigenvalues of the nonlocal IM Hj^. In particular, we find 

?^Mli^n^,^.'-^. (3.20) 

On the other hand an exact expression for i/^"""^-' for all values of p is known explicitly 
from ( p36| ) 

^1 " r(i-/32) r(i-/?2 + 2p)' ^^-^'^ 

where we have used the known vacuum eigenvalue 

M^ac) _ 47r^ r(l - 2/3^) 



of the nonlocal IM Gi [jT| and (|2.36|) . Comparing the asymptotic of ( |3.21j ) with ( |3.20|) one 



obtains exactly the coefficients M which enters the leading asymptotic behaviors ( |2.51J ) 
of the operators A (A): 

M = fr(l - . (3.23) 



Calculating further coefficients in (|3.19| ) we find the large p asymptotic of the eigenvalues 

Tjivac) 

H'n %^+oo {f3 ) P ^. (3.24) 

It is not obvious at all how this asymptotic can be obtained directly from the definition of 
nonlocal IM in terms of the ordered integrals ( |2.16| ) through ( p.l5| ) and (|2.34| ). 



The large p asymptotic of vacuum eigenvalues of the dual nonlocal IM ( |2.37|) is ob- 
tained from (|3.24|) by the substitution P — ^ P 



~( ^, r(n/3-2) r( - i +n(l -/3-2)) , ..-^--i i-2n+^ 

2 y'nnl 

For completeness we also present here the large p asymptotic of the eigenvalues of the local 
IM 

4r-1lp^+oc~(/3^)-^P^-, (3.26) 
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which follow directly from the definition of the operators I2n-i5 see [Q. Note that this 
asymptotic hold for any eigenvalues (not just the vacuum ones) of the local IM. 



Let us now consider the large A behavior of ( 3.18|) . To obtain large A asymptotic 
expansion of A^'"°'^\X) one can close the integral in (|3.18|) in the upper half-plane v > 
—1. Then (|3.18|) can be represented as a sum of infinitely many terms associated with the 



residues of the poles of the integrand located at u > —1, plus the integral over the large 
circle. Note that the contribution of the large circle in fact diverges. Correspondingly, the 
sum of the residues gives only asymptotic series expansion for A('^"=)(A), with zero radius 
of convergence. This asymptotic expansion has the form 



A(""^HA)^C(-)(/3^p) i-y^'^ exp^J^S. {-y')^ 47-1 

ln=0 

{+00 ^ 
-j2 (-1)" {-y^y^ H^r'' 
n=l 

where the variable y is related to A as in ( |2.35| ), 

C'""'' f ) U+„ -i^'y^Ci) (3.28) 

and 

^-- 2jit0^)nS Kl^) ^2^) (^^)^ ■ (^-2^' 

Strictly speaking, here we have derived the expansion ( |3.27| ) in the limit p +oo, so the 
quantities l27-\ Hn""^^ in ( p.27|) denote the p +oo asymptotic of the corresponding 



(3.27) 



vacuum eigenvalues given by ( p.26|) and (|3.25| ). However in the next section we adopt 



additional analyticity assumptions about the operators A (A) which allow for derivation of 
the same asymptotic expansion ( |3.27[ ) for arbitrary p. Moreover, it is natural to expect 
that the expansion (|3.27| ) holds not only for the vacuum eigenvalues but for the whole 
operator A(A); in that case of course the vacuum eigenvalues l27-i Hn"''^^ in ( |3.27| ) 
must be replaced by the operators l2n-i and themselves. 
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4. Conjectures: exact asymptotic expansions and duality 

Consider the following operator-valued function of the complex variable v 



2^ (t{1-(3^) 



^^''^ T{iv^/2) r( - 1/2 + iu/2) r( - + 0/2) (4.1) 
f ^(-A^)-(^+^V^ logA(A). 

J — oo ^ 

The function ^(z/) provides a continuous set of IM parameterized by the variable v. Indeed, 
the operators A (A) with different values of A as well as their arbitrary linear combinations 
commute among themselves and with all local and nonlocal IM. As we shall see below 
the function ^(z/) is remarkable in many respects, in particular, it can be thought as an 
analytic continuation of the local IM \2n-1 to arbitrary complex values of their index n. 
The integral (|4.1D converges only for 2/3^ — 2 < '^mu < —1, however the definition of 



^{ly) can be extended to the whole z^-plane by means of the analytic continuation. The 



latter can be done in various ways. For instance, using the product representation ( p^ ) 
which holds for any eigenvalue A{X) of the operator A(A) one can write the corresponding 
eigenvalue of ^'(z^) as a generalized Dirichlet series of the form@ 



ni^) = 2v^ , / ) , y ( Ai ) . (4.2) 



k=0 



For < < 1/2 this series converges absolutely for 53m u < —1 so that ( [4.2| ) defines an 
analytic function of u in the half plane u < —1 which can be then analytically continued 
to the whole complex z^-plane by the standard technique of the analytic continuation of 
Dirichlet series [^. We expect that the function ^{u) defined in this way enjoy the 



following remarkable analytic properties 



^ In writing (4.1) we make a technical assumption that all zeroes Afe of all eigenvalues of A are 
away from zero, see the remark in Sect. 3. 

^ The definition ( [l.t]) contains an ambiguity in the choice of the branches of the logarithm. In 
(4.2) this translates into the choice of the phases of the roots A^. This ambiguity, however, does 
not afi^ect the following arguments. For definiteness one may assume that | arg A^ | < tt when A^ 
lies outside the negative real axis of A'^, while for real negative A^ the corresponding entry in the 
sum (4.2) is replaced by cosh [ttu{1+^)/2) \Xk\'"''^^+^\ 
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Conjecture 1. For real p and < < ^ the function ^{y) is an entire function 
of the complex variable u. 

These simple analytic properties makes ^(z^) extremely convenient for studying eigen- 
values of the A and T operators. Converting the integral transform in ( f4.1| ) one expresses 
A(A) as 

logA(A) = --^ / — r(l-zz/(l + 0/2)r(zz.e/2)r(-l/2 + zz./2)x 
47r2 Jc^ ^ 

(r(i-/3^))-""%M (-A^)""'«"^ 

where the integration contour is the same as in ( ^.121 ). At the same time, the corre- 
sponding eigenvalue of the transfer matrix T(A) reads 

T(A) = A(qH) + A-^Q-H) , (4.4) 

where 



, w,, -P^ / r(l-zKl+0/2) r(-l/2 + zz./2) 
log A A = 2mP + / ^ ■ cio\ ^ 



(4.5) 



(r(i-/3^))^^™*(.) (_A^)-(^+^)/^ 



The values of the function \I/(z^) at special points on the imaginary z^-axis (where the 
Gamma-functions in ( [4.3| ) and (|4.5|) display poles) are of particular interest. For example, 
it is not difficult to see that the values \1/ ( — 2m (1 — /3^)), n = l,2,...,oo, are related to 
the eigenvalues of the nonlocal IM H^. Using ( p.4|) and (|2.34|) one obtains for the latter 



/ \ —In 

Hn = [r^T{l - /32) j ^ A-2" . (4.6) 
Then it follows from (|4.2| ) that 



fc=0 



-2n 



For other special values on the imaginary z^-axis we will adopt the following 

Conjecture 2. For realp and < < | the operator ^{y) has the following special 
values on the imaginary u-axis 

*((2n-l)z) = l2n-i , n = 0,l,..., (4.8) 
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-+1 



*(2m(/3-2 -1)) 



*(0) = P , 



(4.9) 
(4.10) 



where l2n-i o.nd H„ denote the local IM (\2. 22^ ) and the dual nonlocal IM (\2.3'/[ ) respec- 
tively. 



Note that (|4.8[ ) with n = reads 



*(-z) = I_i = I 



(4.11) 



where I is the identity operator. The relations (|4.8|) - ([4.10D together with the conditions 
( p.23|) , (|2.37|) makes it natural to assume that the operator ^ satisfy the following duality 
condition: 

Conjecture 3. 



(4.12) 



An initial motivation for the above conjectures came from the study of the large p 
asymptotic of the vacuum eigenvalues of A(A) in the previous section. Indeed, comparing 
( p.l8|) and ( [1.3|) one obtains 



(4.13) 



which with an account of ( |3.26|) and (|3.25| ) obviously satisfy all the statements of the 
Conjectures 1-3 i. Further motivations and justifications of the these conjectures are 
discussed below. 

We are now ready to derive exact asymptotic expansions of the eigenvalues for large 
A^. This is achieved by calculating the integrals ( [4.3| ) and ( |4.5D as formal sums over residues 
in the upper half plane Qmu > — 1. Using (|4.8|) one thus has from ( [1.5|) for — » oo 



log A(A) ~ z m I ( - A^ ) 
where 



2-2/3^ 



1 — 2n 

5^(-l)" a (-A^)^^ ^(l-2n)(l+0 ^ (4.14) 



n=l 



m 



2v^r(i-|) 
r(i-|) 



(r(i-/3^ 



1+^ 



(4.15) 



5 

ones. 



Like in (3.26) , the asymptotic ( [4.13| ) holds for ah eigenvalues of not just the vacuum 
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and 



r(i + (n-i)0 



(4.16) 



It follows then that the large A asymptotic behavior of the operator T is given by ( |2.21| ) , 
where the numerical coefficients m and C„ are given by the formulas (|4.15 ) and ( [4.16|) . 
The asymptotic expansion (|2.21|) holds for 



— 7r < arg A < tt . 

Similarly one can get the asymptotic expansion for the operator A(A) 
operator 

2 P 



-2P 



(4.17) 



Define the 



(4.18) 



which does not depend on the spectral parameter A. It is also convenient to use the variable 
y given by ( |2.35|) instead of the variable A and exhibit all arguments of A{iy): 



A{zy,/32,(^(«)}^exp(-5^ (-1)- (y')"H, 
Then, from (|4.3|) one has for \y\ 00, —tt < arg y'^ < % 

{00 
Y^B^iy^)^ hn-i 
n=0 

A{zy-i^,r',/3-VH} , 
where the coefficients are given by ( p.29| ) and 

00 

A{zy-i^,/3-',rVM}=exp(-5^ (-1)- [y^yP H 



X 



(4.19) 



(4.20) 



n=l 



The explicit form of the operator C in ( |4.19| ) is not determined by the above calculations. 
For the vacuum eigenvalue its p —>■ +00 asymptotic is given by (|3.28|). Moreover, one can 



show that [24 



C(""'^)(/32,0) = . (4.21) 

We have the following conjecture about the exact form of this coefficient for all values of 

1 
2 

, . . , — V(^-\-2r)B-^ 

(4.22) 



p and < < 1 



r(i + 2p) 
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The formula for the asymptotic expansion (|4.19 ) is essentiaUy equivalent to the Con- 
jectures 1,2 given above. In fact, the coefficients in front of powers of A in ( [4.19| ) are in 
one-to-one correspondence with the values ( [4.8| ) while any singularity of ^(z^) in the upper 
half plane would bring in some additional terms in ( [4.19| ). 

The asymptotic expansions (|2.21|) and ( |4.19|) are in remarkable agreement with the 
numerical calculation through the (modified) TBA equation. We postpone a detailed 
description of these calculations to a separate publication but just mention some the results 
here. In we calculated numerically a few coefficients in (|2.21| ) for two vacuum states 
(A = —1/5 and A = 0) in the Ai2,5 CFT (c = —22/5) and found an excellent agreement 
with the corresponding exact eigenvalues of the local IM given explicitly in [|I| (up to 
I7 inclusive). Another numerical result concerning the part of ( [4.19|) containing the dual 
nonlocal IM is mentioned in Sect. 5. 

As an additional support to our conjectures consider the the case 

2 ' 

where the eigenvalues of T and A-operators can be calculated explicitly. Moreover, the 
eigenvalues of the local IM can be independently found using the fermionic representation. 
The value = | does not lie in SD (|2.18| ) therefore the results of Sect. 2 do not apply 
directly. In particular, the definitions ( |2.11J ) and ( |2.29|) for T(A) and A(A) requires a 
renormalization since the nonlocal IM in ( |2.15| ) and (|2.33| ) diverge logarithmically at 
= i. It turns out that this renormalization affects the functional equation ( |3.3| ). To 
see this consider the expressions for the vacuum eigenvalues to within the first order in A^ 

j.{vac) ^ 2 cos(27rp) + A^ gS""'^ + 0(A4) , (4.23) 

A^-^HX) = 1 - f ( g\-^^ - sin(27rp) ) + ^(A^) , (4.24) 
4cos 27r» V / 



where 



gS""''^ = 2 cos(27rp) ( 2C - ttV' (- - 2p) - TTi/; (- - 2p) ) , (4.25) 

il}{x) = dj;\ogT{x) is the logarithmic derivative of the gamma-function and C is a (non- 
universal) renormalization constant depending on the ultraviolet cutoff. A simplest way 
to obtain these expressions is to set = | — e, e — in (|3.21| ) and (|3.22| ) ; that gives an 
analytic regularization of the divergent integrals with the value of C 

C=^ + 7r^(l). (4.26) 
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With the above accuracy in the eigenvalues ( f4.23| ) (|4.24 ) satisfy a "renormahzed" func- 
tional equation 

T{X)A{X) = e'^'^'P-'^^^^ A{q\) + e-2"^P+*'^'^'A(g-U) , (4.27) 

where q = exp(i7r/2). Using the lattice regularization (i.e., considering discrete approxima- 
tions to the P-exponents in and ( |2.29| ) and then tending the number of partitions 



to infinity) one can show that the functional equation ( [4.27| ) is, in fact, exact in the sense 
that it is valid for arbitrary eigenvalues T(A) and A{X) to all order in A^. 

The functional equation (|4.27|) completely determine the eigenvalues T(A) and A{X) 



provided one assumes them to be entire functions of A^. One can obtain 

L 

T(A) =T(^-)(A) l[F{iX,p,nt,n]:) F{i\ -p,nl ,nt) , (4.28) 



fe=i 



A(A) = )(A) n^(A,P,n+,n-) , (4.29) 



where 



k=i 



+ _ (2p - n- + I - 7rA^)(2p + n+ - |) 



F(A,p,n+,n-) = ^ —-j^^ (4.30) 

and n^, n^, ■ ■ • , are two finite sequences of non-negative integers, 1 < nf < ■ ■ ■ < n^, 
which uniquely specifies certain vector in the Fock space J^^; this vector has the Virasoro 
weight 

L 

A{p; nt,...,n+-n^,...,nl} = A^"-) (p) + 5^(n+ + - 1) . (4.31) 

k=l 

For the vacuum eigenvalues L = and 

^ ' r(l/2 + 2p + 7rA2) r(l/2 -2p + ttA^) ' ^ 

A(-) (A) = e-''' ^i^P+y^) ^4.33) 
^ ^ r(2p-M/2-7rA2) ^ ^ 



^ The basis of this kind in the bosonic Fock space is usually referred to as "fermionic basis" 
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It follows then from (14. 21) that 



(4.34) 



where ({s, a) is the generalized zeta function defined as the analytic continuation of the 
series 

oo 

C(s,a) = ^(a + n)-" , 5?e s > 1 (4.35) 

n=0 

to the whole complex plane of the variable s. This function is analytic everywhere in the 
s-plane except the point s = 1, where it has a simple pole. Therefore the eigenvalues ( [4.34| ) 
are entire functions of u in agreement with our conjecture. 

The values of ( [4.34| ) at the integer points on the imaginary z^-axis follow from the 
formula 

Bm+i{a) 



C(-m, a) 



m + 1 



m = 0, 1, . . . 



where Bmio:) are the Bernoulli polynomials [^. Note, in particular, that 

^'(-z) = l, ^{0)=p, 



(4.36) 



(4.37) 



in agreement with ( [4.11|) and (|4.1C1| ). Further, the values of ( [4.34|) at u = z(2n — 1), 
71 = 0, 1 . . ., conjectured in ( [4.8| ), perfectly match the eigenvalues of the local IM, which 
for = ^ can be independently obtained from the explicit expression for the local IM 
through free fermion fields 



j{vac) 
hn-1 



2-" B2n{2p+^) . 



(4.38) 



The equation (O) can not be tested since the eigenvalues of the dual nonlocal IM are 



not generally known independently except for the vacuum eigenvalue of H\ 
latter one we have 

vl>M(2z) = 2-3^3(2^+1) =p (p-i) (p + i) . 



iyac) 



For the 



(4.39) 



in a precise agreement with ( [4.9| ) and the relation which is dual to ( |3.21| ). 

Finally note that the asymptotic expansions of ( [4.28| ) and ( ^.29| ) agree with ( |2.21| ) 
and ( [4.19| ) provided one calculate the latter in the limit = ^ — e, e ^ 0, and identifies 
the divergent part of the coefficient in the leading asymptotic of ( p.21|) and (|4.19|) with the 
renormalization constant C in ( [4.26| ), (|4.32| ) and (|4.33| ). 
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5. Vacuum eigenvalues of the Q - operators and non-equilibrium states in 
boundary sine-Gordon model 

As is explained in Sect. 2 the vacuum eigenvalue A*^""^^(A ) for integer 2p = N coincides 
with the Coulomb gas partition function ( |2.5UD . In fact it was already observed in 
that the Coulomb gas partition function ( p. 501 ) with p = satisfies the functional equation 
( p..4|) B. The Coulomb gas partition function is obviously related to the finite temperature 
theory of mass-less bose field x) on the half plane x < with interaction at the 
boundary; its action is 

^=^j_ dtj_ rfx($2_$2)^_ dtcos ($(t,0) + W) . (5.1) 

Here k, and V are parameters i. This model finds interesting applications in dissipative 
quantum mechanics [Q, As was discussed in [0, [|4| it also describes 

the universal current through the point contact in quantum Hall system. At nonzero 
driving potential V and arbitrary temperature T the system ( |5.1D develops a stationary 
non-equilibrium state with 

Jb = ^{^.{t,x)) = -K {sin {^{t,0) + Vt)) ^ 0. (5.2) 

This quantity is interpreted as the backscattering current through the point contact i. 
In general case (|5.2|) is non-equilibrium expectation value and as such it requires non- 
equilibrium methods for its computation |T^. Note that if the driving potential is continued 
to pure imaginary values 

V = Vn = 2niN p-'^T (5.3) 



^ We obtained this result independently, along with more general statement ( p. 31 ), before the 
paper [32| appeared. 

^ The composite field cos{^ + Vt) in ( p.lD is assumed to be canonically normalized with respect 
to its short-distance behavior, i.e. cos ($(t) + Vt) cos ($(t') + Vt') ~ 2"^ (z(t - t') + )"^'^ , as 
is conventional in conformal perturbation theory. This is why no ultraviolet cutoff will appear in 
the matrix elements below. 

^ The voltage V and current Jb ( ^.2[) differs in normalization from the real voltage V^^'^^'^^ 
and total current j'^^^y'^^ in the Hall system 

y(phvs) = g-l y ^ jiphys) = ^ jf {v + Jb) , 

h ^ ' 

where e and h are the electron charge and Plank's constant. Also, 0^ coincides with the fractional 
filling of the Luttinger state in a Hall bar and the temperature is measured in energy units. 
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with integer A^, one can make the Wick rotation t — > —ir in ( ^.11) and formaUy define the 
Matsubara partition function 



exp(-^M) , (5.4) 



where 

^M = -r^l dr f dx{^l + ^l)-^[ dr cos {^0,t) + 27VN (3-^Tt) 

4'iTp Jo J-oo P Jo 

(5.5) 

and the functional integral is taken over Euclidean fields $ which satisfy the Matsubara 
condition $(r + (3'^/T,x) = $(r, x). It is easy to see that up to overall constant (the 
partition function (|5.4|) with k = 0) ( ^.4|) coincides with the series ( p.50|) with 



2T V 27rT 

Moreover, it is possible to show, using non-equilibrium methods, that for V = 
the expectation values of e^***-*'^) in (|5.2| ) can be calculated in the Matsubara theory 
(U) as 

(e±** ) = J e±'*("'°) e-^^' / Zn{i^) ■ (5.7) 

Of course equilibrium state can not support nonzero current and indeed it is not difficult 
to show that 

{e^-(e^'^)-e-^-(e-^^)}| =0. (5.8) 



Nonetheless it is natural to expect that the non-equilibrium expectation value ( |5.2|) can 
be obtained by some kind of analytic continuation of ( |5.7| ) back to real V . Needless to 
say this analytic continuation is ambiguous. 

In |T^ the current ( |5.2| ) is calculated exactly (for integer values of /3~^) using the 
Boltzmann equation with the distribution function of charge carriers determined through 
Thermodynamic Bethe Ansatz technique. A conjecture is proposed about exact current 
for arbitrary and V in the recent paper [|TB[. It was suggested there to define the 
"partition function" Z2p{p) for V = Ainp /?~^T as certain analytic continuation of 
which uses infinite sum expressions for the coefficients in ( p.50| ) obtained with the help of 
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Jack polynomials It is conjectured in that the current ( p.2|) can be expressed 

in terms of Z2p{n) as 

JB{V,^,(3')=^nT^^,logp^^, V = ■ (5.9) 



This conjecture agrees with the earlier conjecture for the linear conductance in In 
fact, various checks performed in fl^ suggest that Zip thus defined satisfies the functional 
relation (|3.3|) . It is therefore very plausible that Z^p coincides with the vacuum eigenvalue 

Indeed, we found that the result of is in complete agreement with the following 
formula 

Jb(F, ^, = mT A^Alog + . , (5.10) 

A™ (A) 

where 

sin(7r/32) z/32y 
TT 47rJ 



Namely, the expression for Jb obtained in [|T5| can be written in the form ( 5.10| ) 



with A±{X) solving the functional equation (p.3|). Note that for 2p = N both A^]^"''^^ and 
^{"vo-c) (,Qj]^(,j^g with Ztv, so each of them defines certain analytic continuation of Z^r to 
complex A^. It is suggestive to note though that analytic properties of A^"'^'' and 
as the functions of p are different, namely is analytic at 5Re 2p > — whereas 

aL""''^ is analytic at 3?e 2p< (3^. There is no clear notion of partition function for a non- 
equilibrium state, and therefore it is remarkable that A^^"'^^ (and A^_^"'^^ as well) admits 
interpretation as an "equilibrium-state" partition function of the system similar to 
but with additional boundary degree of freedom described by the "q-oscillator" £±,7i in 

(Pi). 

According to ( [4.19| ) the conjecture (|5.10| ) implies in particular that the backscattering 
current (5^) satisfies the following "strong- week barrier" duality relation 



Jb{V,ix,(3'')^-V -[5-'' JB{f3''V,Cii ^, (5.11) 
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We acknowledge a private communication with P. Fendley, F. Lesage and H. Saleur who 



explained to us how the definition of Z2p in [16| must be understood. 

Again, we have arrived at (5.10) independently, before the paper [16| appeared. However 



we were significantly influenced by the conjecture about the universal conductance proposed in 



[24] and by the results of [15|. 
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(remarkably, the factor containing the local IM in ( 4.19 ) cancels in the ratio 

which generalizes similar relation obtained in []TB| for the case The constant 



C = C{/3^) in ([5TTTD reads exphcitly: 



(5.12) 



Conventionally, one introduces the renormalized coupling parameter X |]T^, p| 



^2 .1-2/5^ /- r(i + /3^) 2 



(5.13) 



defined in such a way that 



dvJB{V,X,/3'')\y^^ = -X^ + 0{X^) , X' 



(5.14) 



According to the formula ( p.llD the current admit the following decompositions: 

oo 

Jb{V,X, {3')=-J2Fn{V,f3')X'- 

n=l 

oo 



(5.15) 



2n 



where are certain functions and the constant K = K{I3'^) is 



vr\^ + ii^ r r3(l + /32) r3(l+/3-2) 



r(i/2 + /32) J r(i/2 + /3-2) ■ 



(5.16) 



Note that this constant determines the leading X ^ oo asymptotic of the backscattering 
current, 

JB^-V + 2r'TK{P^) sinh^ r(r^ + z^^)rjr^-z^) ^-^^^^^-^^ _ 

(5.17) 
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In fact, exactly this duality relation (for the nonlinear mobility in associated dissipative 



quantum mechanics problem) was proposed a while ago in [35|. The arguments in |35| are based 
on "instanton" description of the hopping amplitudes in the strong barrier limit; from general 
point of view this description could be regarded as an approximation. Therefore it looks quite 



remarkable to us that this relation is indeed exact. See also the discussion in [15|. 
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We should stress that the formula ( 5.10 ) is a conjecture and therefore any checks 
would be valuable. As was mentioned above this formula agrees for integer f3~'^ with the 



TBA solution obtained in ||T^. Although the calculations in are also based on a 
conjecture (exact applicability of the Boltzmann approximation to integrable systems), 
one could check ( p.lO|) against numerical results obtained from the TBA solution. We 
have done that for = 3 and found an excellent agreement with ( |5.17D . In particular 
our numerical value for the constant K is 

i^n^m(l/3) = 3.35485280612... , (5.18) 

(note that this number is slightly off from the numerical result for the same constant given 
in 1^^) which is in agreement with the exact value 

K{l/3) = bl^Ll l^^li^^^ ' = 3.35485280611990... . (5.19) 

Completely rigorous check can be made in the case = 2, where the formula ( |5.1Q ) 
agrees with the explicit calculations in the free-fermion theory [|l^, [|15| 

^-(^'^'2) = -^(^(2 + ^ + 4^)-^(2+^-4^))' ^'-''^ 

where 'il^{x) = dxlogV^x). Another interesting limiting case is the classical limit /32 — > 0. 
In this limit the eigenvalues A^^'^'^\x) reduce to (see Appendix B): 



A^r^\X)\^.^o ^ r(l ± p) J±p{V2X) , (5.21) 

where the variables 



X = 72 /3-2 A , p = 2/3-2 p 

are kept fixed when (3'^ 0. Substituting this into ( |5.10| ) one obtains after a little algebra 

/ X 2Tsinh77^ iV 
Jb(V, X, 0) = -V + ^ — p = 5.22 



where Ip{x) is modified Bessel function. As is known the field theory ( |5.1J ) can be inter- 
preted in terms of dissipative quantum mechanics of a single particle in periodic potential 
g, 0, [|5[, [||. As /32 plays the role of the Planck constant, in the hmit ^ this 
reduces to the theory of classical Brownian particle at finite temperature T. In Appendix 
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C we study associated Fokker-Planck equation and calculate the classical current. This 



gives additional support to the conjecture (|5.10|) . 

As was explained above, the vacuum eig envalue A(""^)(A) at p = coincides with the 
partition function of (|5.1| ) at the temperature T and zero voltage V = 0. Therefore the 
asymptotic expansion ( |4.19D , ( ^4.2U| ) specialized for the vacuum eigenvalue with p = is 
essentially the low temperature expansion for the associated free energy. In particular, the 
leading asymptotic of the heat capacitance C ( X, ) of the point contact in the quantum 



Hall system at T ~ const X ^-p'^ and zero voltage V reads 



v/^/34r(2-^) I 2r(i+/32) / ' 3 

(5.23) 

Likewise, ( |2.21|) determines the low temperature expansion for the impurity free energy in 
the s = 1/2 anisotropic Kondo problem. 



6. Discussion 

In this paper we have studied further how the powerful apparatus of the Yang-Baxter 
theory of integrability can be brought about directly in continuous Quantum Fields Theory. 
We have constructed the operators Q±(A), which are field theoretic versions of the Q - 
matrix of Baxter, directly in Conformal Field Theory. The Q - operators are constructed 
as the traces of certain monodromy matrices associated with the infinite - dimensional 
spaces - the representation spaces of "q-oscillator algebra" (|2.27|) . It is worth mentioning 
that our construction is not specific for the continuous theory - the Q- matrix of lattice 
theory admits similar representation (we will give the details elsewhere). We also found 
that the Q - operators thus constructed satisfy the remarkable relations ( pj.40| ), ( |2.41| ), 
( p.42|) and ( |2.43|) . These relations allow one to employ the powerful machinery of nonlinear 
integral DDV equations to study the eigenvalues of the Q - operators in the highest weight 
Virasoro module. We have used the DDV approach to derive (under some analyticity 
conjectures) various asymptotic expansions for both Q and T operators. We also observed 
a remarkable (although somewhat puzzling) relation between the vacuum eigenvalues of the 
Q - operators and the stationary transport characteristics in boundary sine-Gordon model 
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(the later also relate to the kinetic properties of one-dimensional quantum particle coupled 
to a dissipative environment) . In this paper we did not present the derivation of our basic 
relations (|2.40| ), ( |2.41| ), ( |2.42|) and (|2.43|) ; this gap will be filled in the forthcoming paper 

in. 

Clearly, further study of the Q and T operators is desirable. First, almost all the 
discussion in this paper concerns the "Semi-classical Domain" < 1/2 (the case = 1/2 
is studied explicitly through the free fermion theory in the Sect. 4). However, the most 
interesting CFT (notably, the unitary CFT) lay outside this domain. We have argued in 
the Sect. 2 that the Q operators can be defined for 1 > > 1/2 as well, but there are 
reasons to believe that outside the SD the analytic properties of both Q and T operators 
undergo significant (and very interesting) changes. We are planning to extend our analysis 
to the domain 1 > > 1/2 in the future. 

The analysis in this paper is concerned explicitly with the conformal field theories 
(more precisely, the chiral sectors of a CFT). One can argue however that the operators 
Q-I-(A) (as well as the T operators of |^]) can be defined in non-conformal integrable QFT 
obtained by perturbing the CFT with the local operator \I/i,3, the most important relations 
( p.40|) , (|2.41| ), ( p.42|) and ( p.43|) remaining intact, and the most significant modifications 
being in analytic properties of these operators (they develop essential singularities at 
in the perturbed case). The Q operators and associated (modified) DDV equations can 
be used then to study the finite-size spectra of these non-conformal theories. We will report 



some preliminary results in this direction in the forthcoming paper ||3q] . 

And finally it seems extremely desirable to get more understanding about the relation 
of the Q - operators to the non-equilibrium properties of the boundary sine-Gordon model 
discussed in Sect. 5 above. In particular it seems important to find a physical interpre- 
tation to the "g - oscillator" degrees of freedom which evidently play central role in our 
construction of the Q - operators. 
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7. Appendix A. 



In this Appendix we give details of the derivation of the DDV equation (|3.8|) . For 



any eigenvalue A{X) satisfying (|3.3|) consider the function a(A) defined by (|3.6|) . For real 
values of p the property (ii) of the function A{X) (given in Sect. 3 above the formula (|3.4[) ) 
implies 

a{Xr=a{X*)-' , (7.1) 

where the star denotes the complex conjugation. Next, it follows from leading asymptotic 
( CT ) of A{X) at large A that 

loga(A) ~ -2iM cos(7r^/2) A^+\ A^ ^ oo , | arg A^ | < 27r/32 , (7.2) 

where = Moreover, the function a(A) obviously remains finite for small A^ 

a(A) -47rzp + 0(A2), A^ ^ . (7.3) 

The function a (A) satisfies the Bethe-Ansatz type equations (p.7|). The product represen- 
tation (^) implies 

oo 

loga(A) - imp = F{XX-^) , (7.4) 

where 

nX) = log ■ (7.5) 
The sum in the RHS of (|7.4D can be written as a contour integral 

loga(A) - Amp = /(A) + [ ^ F{X^i-') c»>g(l + a(//)) , (7.6) 

where 

/(^) = T!nXX-') (7.7) 

a 

denotes a finite sum including only those zeroes A^ which do not lie on the positive real 
axis in the A^-plane. The contour C goes from +00 to zero above the positive real axis, 
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then winds around zero and returns to infinity below the positive real axis in the A^-plane. 
Integrating by parts (boundary terms vanish due to ([7 .21), since < ^ < 1 for in SD 
( p.lSp ) one obtains 

0=/(A)- r Xd^F{X^i-') {log(H-a(^ + zO))-log(l + a(^-zO))} 

= / ^ \dxF{\^x-^) 53mlog(l + a(;u-z0)) - / \dxF{\^-^) loga(;u) 



(7.8) 



where we have used (Ol)- 

Introducing new variables 6*, Q' and as 



A = ei+«, ;U = ei+«, Afc = ei+« (7.9) 



ITT ^ 

and recalling that q = 6^+^, one obtains 



de' R{6 - e') 9m log (1 + a{e' - zo)) + / de' R{e - e') log a{e') . 

-oo J — oo 

(7.10) 



where 



With the standard notation for the convolution ( p.lOj) , the equation ( [7.10|) can be written 



as 

K*\oga{e) ^ Anip + f(e) -2i R * Qmlog {l + a{6 ~ iO)) , (7.12) 

where 

K{9) =6{9) - R{9) . (7.13) 



Let us now apply the inverse of the integral operator K to both sides of ( [7.12 ); in this one 



has to add an appropriate zero mode of K to the r.h.s. of the resulting equation to make 
it consistent with the asymptotic conditions ( |7.2| ) and ( [7.3) ). In this way one obtains the 
integral equation (|3.8| ) in the main text. As is noted there, this integral equation has to 
be complemented by a finite number of the transcendental equations ( p77| ) for the roots A^ 
lying outside the positive real axis of A^. 

For the vacuum eigenvalues with 2p > — all the roots A| are real and positive and 
the term i logS{6 — 6a) in the r.h.s. of ( p.8|) is absent. Let Ag be the minimal of the 
zeroes A^. Introduce the function a(""^)(A) related to A(""'^)(A) as in ( p.6| ) and denote 

B{p)=^-^ logA^ (7.14) 
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The function log (l + a^'"^'^\6)) is analytic in the strip — tt^ < 6 < with the branch 
cut along the positive real axis from B{p) to infinity. So it is analytic for real 9 < B{P), 
where it obeys the relation 

log 0(6*) = 2Qm log (1 + a{6 - iO)) , (7.15) 

which follows from Note that the infinitesimal shift — zO here is not essential. Taking 

( [7.15|) into account one can rewrite the ( p.8|) for 9 < B{p) as 



-| + M cos^ - — de\ogS{e - 9') log (l + a(""^)(^')) 

I'+oo in/ 

= / aelog^(6' - 9') 53m log (l + a(""'^V6'' - iO)) . 

Jb(p) 27r« 



(7.16) 



Let us make technical assumption 



B{p) ~ const log p as p ^ +00 . (7-17) 



Then one can show that the r.h.s. of (|7.16|) decreases at the large positive p and therefore 
can be dropped in the leading approximation at p ^ +00. This brings ( [7. 16] ) to the linear 
integral equation ( |3.14| ) of the Winer-Hopf type. 



8. Appendix B 

In this appendix we consider the functional equations for the eigenvalues of T and Q 
operators in the classical limit 0, where they reveal a remarkable connection with 

the theory of classical Liouville equation 111. We start with the functional equations for the 
eigenvalues Tj{X) of the operators ( ^.11| ) which follow from ( f2.13D ; they can be written in 
the form 

T,{q-^X)T,{q-^X) = l + T^_.{X)T^^^iX) (8.1) 
Consider limiting values of the eigenvalues Tj{X) and Q±{X), when 

2A = /32e^ J = 42' 2p = /32p, /32^0, (8.2) 
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The fact that the functional relations associated with the TBA equations have many features 



in common with the classical Liouville equation is well known to experts |26]. 
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and the variables o", r and p are kept fixed. We assume that in this hmit 



-</)(cr,T) 



(8.3) 



where (j){a, r) and Q±{(j) are smooth functions of their arguments. Then it is easy to see 
that the functional equation (|8.1j ) becomes the classical Liuoville equation 



(a^+a^) <^(a,r)=e 



for the field (j){a, r) in the Euclidean space. It should be complemented by the periodical 
boundary conditions 

(^((T + ZTT, r) = 0(cr, r) , (8.5) 



since Tj{X) in (|8.3| ) is a single- valued function of A^. The limiting form of the functional 
equation ( p.43|) reads 



e-^(a,r) _ (4^p)-i ( Q+(a + zr)Q_(a-zr) - Q+(a - zr)Q_(a + zr) ) 
whereas ( |2.42| ) becomes the ordinary Wronskian condition 



.6) 



(8.7) 



for the function Q±{a). 

One can recognize in ( |8.6|) , (|8.7|) the general local solution of the Liouville equation 
satisfying the condition 



-0(cr,T) 



r + 0(r2) , T ^ 0. 



The Baxter's relation (|1.4|) reduces to the second order differential equation 



with a periodic potential 



w{a + m) — w{a) 



(8.9) 



(8.10) 



which is determined by the leading asymptotic of the eigenvalue T{X) = Ti (A) in the limit 
(H) 

T{X)=2 + {Tr(3f w{a) +0{(3^) (8.11) 
37 



The functions Q±{u) are just two linearly independent Bloch-wave solutions 



Q±{a + in) = e^'^'P Q±{a) . (8.12) 



to the second order differential equation ( |8?9| ) . Note that ( |8.12|) imply the periodicity 
of the solutions ( p^ ) for the Liuoville equation ( |8.4|) . 

It is illustrative to apply the above limiting procedure to the vacuum eigenvalues of T 
and Q operators. Notice that all the vacuum eigenvalues of the nonlocal IM ( |2.16|) entering 
the series expansion ( |2.15D remain finite in the limit ( |8.2|) . Therefore only two first terms 



of the expansion 

j.(.ac) _ 2 cos(27rp) + A^g!"""^ + 0{X^) . (8.13) 



contribute to ( |8.11| ). Using ( p.22|) one thus obtains 

w{a) = e^" - p'^ (8.14) 
The solutions to the differential equation 

(dl + e^'^-p') Q±(a)=0, 
satisfying the conditions ( ^.7| ), ( fj.l2D have the forms 



Q±(p) = r(l±p) J±,(e^) , 

where Ja{x) is the conventional Bessel function. 

The eigenvalues Tj^"'^\x) in the limit (|8.2[) are then expressed from ( |8.3| ) through the 
corresponding solution ( |8.6|) of the Liouville equation 



^ ^ > 71/32 



1 



2i(3^ sinvrp 



.15) 



It would be interesting to check this formula against the solutions to the TBA equations 
in the limit — > 0. 
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9. Appendix C 

As is discussed in the field theory describes a dissipative quantum mechan- 
ics of one dimensional particle if one interprets the boundary value = $(t,0) as the 
position of the particle, the bulk degrees of freedom $(t,a;), x playing the role of 
thermostat. As enters ( |5.1D as the Planck's constant, in the limit — > this prob- 
lem reduces to the one of classical Brownian particle which is described by the Langevin 
equation 

^ m = sin {m + Vt) + ^(t), (9.1) 
where $ = 9t$(t), and ^(t) is the white noise 

{mm) = -s{t-t'). (9.2) 

TT 

We are interested in the limiting value of the average velocity 

JBiV) = {^)t^^ . (9.3) 
The time dependence of the potential term in ( |9.1|) can be eliminated by the substitution 

z{t) = $(t) + Vt , (9.4) 

which transforms ( p.l| ) to the form 

z = V -2nK sm{z) + 2n^ . (9.5) 



The driving force V appears here because the transformation (|9.4|) brings us to the frame 



which moves with the velocity —V with respect to the thermostat. The current (|9.3|) then 



IS 



JBiV) = -V + JiV) , J{V) = {z)t^^. (9.6) 



The probability distribution P{z) associated with the stochastic process ( p.5| ) satisfies the 
Fokker-Planck equation 

dtP = 2tiT d,{{-v + sinz) P + d^P] , (9.7) 



where 



(9.8) 



The suitable solution of ( p.7|) , which describes the stationary drift has the form (see e.g. 

m) 



Jz 27r 



z+27r 

cos(y)-uy 



where 

-2tt 



(9.9) 



J\f= dz Pq{z) (9.10) 
Jo 

is the normalization constant. The current J can be then expressed as 

J = 2nT Af-^ (1 - . (9.11) 

The integral ( p.lO| ) is calculated explicitly with the result 

2T sinh( tti/) , , 

J= ^ ^ — , 9.12 

h.iV2X) I-^,{^/2X) 



where Ip{x) is the modified Bessel function. With the account ( |9.6D , (^.8]), (p?B|), ( [5.13| ) 
this agrees perfectly with ( p.22| ). 
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